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rvQ ' Abstract. We study the boundary symmetries of models of nonequilibrium physics where 

the steady state behaviour strongly depends on the boundary rates. Within the matrix 

product state approach to many-body systems the physics is described in terms of matrices 

defining a noncommutative space with a quantum group symmetry. Boundary processes 

OO ' lead to a reduction of the bulk symmetry. We argue that the boundary operators of an 

interacting system with simple exclusion generate a tridiagonal algebra whose irreducible 
representations are expressed in terms of the Askey- Wilson polynomials. We show that the 

i-C . boundary algebras of the symmetric and the totally asymmetric processes are the proper 

limits of the partially asymmetric ones. In all three type of processes the tridiagonal algebra 
r^ . arises as a symmetry of the boundary problem and allows for the exact solvability of the 

model. 
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SD. ' 1 Introduction 

'NT 



Out of the rich variety of phenomena in nature the most interesting occur in nonequilibrium 



o 

QQ ' conditions and their complex behaviour is far from being well understood. The study of nonequi- 

^D ■ librium phenomena is of current interest and the way to describe the general characteristics of 

a system out of equilibrium goes through the analysis of mathematical models. Such models 
have to be simple enough, still they must be physically significant to exhibit the structure of 
the complex phenomena. Stochastic interacting particle systems [H El [3] received a lot of atten- 
tion. Among these, the asymmetric simple exclusion process (ASEP) has become a paradigm 
in nonequilibrium physics due to its simplicity, rich behaviour and wide range of applicability. 
Introduced originally as a simplified model of one dimensional transport for phenomena like 
kinetics of biopolymerization [6], it has found applications from traffic flow [H [5]) to interface 
growth [7], shock formation, hydrodynamic systems obeying the noisy Burger equation, problems 
of sequence alignment in biology [8] . 

At large time the ASEP exhibits relaxation to a steady state, and even after the relaxation it 
has a nonvanishing current. An intriguing feature is the occurrence of boundary induced phase 
transitions [2] and the fact that the bulk properties depend strongly on the boundary rates. 

The asymmetric exclusion process is an exactly solvable model of a lattice diffusion system 
of particles interacting with a hard core exclusion, i.e. the lattice site can be either empty 
or occupied by a particle. As a stochastic process it is described in terms of a probability 
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2 B. Aneva 

distribution P{sm, t) of a stochastic variable Sm = 0, 1 at a site m = 1, 2, . . . , L of a linear chain. 
A state on the lattice at a time t is determined by the set of occupation numbers si,S2, ■ ■ ■ ,sl 
and a transition to another configuration s' during an infinitesimal time step dt is given by 
the probability r(s,s')dt. With the restriction of dynamics that changes of configuration can 
only occur at two adjacent sites, the rates for such changes depend on these sites only. The 
two-site rates T(s, s') = p'^^'^'^+i = Yl^j, i,j, k,l = 0,1 are assumed to be independent from the 

position in the bulk. Due to probability conservation r(s,s) = — X^^'-^s r(s', s). For diffusion 
processes the transition rate matrix becomes simply F^^ = gi^- At the boundaries, i.e. sites 1 
and L additional processes can take place with rates L^ and R^ {i,j = 0,1). In the set of 
occupation numbers {si,S2, ■ ■ ■ , sl) specifying a configuration of the system Sj = if a site i is 
empty, Sj = 1 if there is a particle at a site i. Particles hop to the left with probability goidt 
and to the right with probability giodt. The event of exchange happens if out of two adjacent 
sites one is a vacancy and the other is occupied by a particle. The symmetric simple exclusion 
process is known as the lattice gas model of particle hopping between nearest-neighbour sites 
with a constant rate ^fQi = gio = g. The partially asymmetric simple exclusion process with 
hopping in a preferred direction is the driven diffusive lattice gas of particles moving under 
the action of an external field. The process is totally asymmetric if all jumps occur in one 
direction only, and partially asymmetric if there is a different non-zero probability of both left 
and right hopping. The number of particles in the bulk is conserved and this is the case of 
periodic boundary conditions. In the case of open systems, the lattice gas is coupled to external 
reservoirs of particles of fixed density. The most interesting examples (see [9] and references 
therein for a review) are phase transitions inducing boundary processes [H [10] when a particle 
is added with probability adt and/or removed with probability ^dt at the left end of the chain, 
and it is removed with probability (3dt and/or added with probability 5dt at the right end of 
the chain. 

The time evolution of the model is governed by the master equation for the probability distri- 
bution of the stochastic system — ^^ = ^^/ F(s, s')P{s' , t). It can be mapped to a Schrodinger 
equation in imaginary time for a quantum Hamiltonian with nearest-neighbour interaction in 
the bulk and single-site boundary terms Jj. ' = —HP(t), where H = ^ • Hjj^i -\- H^^' + H^^> 
and the ground state of the Hamiltonian, in general non-Hermitian, corresponds to the steady 
state of the stochastic dynamics where all probabilities are stationary. The mapping F = 
-VqU'^HxxzU^, where 

(see jll] and [12] for the details) provides a connection to the integrable SJ7q(2)-symmetric XX Z 
quantum spin chain with q = ^; Hxxz is the Hamiltonian of the Uq{su{2)) invariant quantum 

spin chain Hxxz 1^3] with anisotropy Aq and with added non diagonal boundary terms Bi 
and Bl 

Hxxz = H-^j^^ + Bi + Bl 

L-\ 

= -1/2 Y^iafaf^, + afaf^, - A,af<+i + h{a^^, - af) + A,) +B, + B^. 

The transition rates of the ASEP are related to the parameters A^ and /i, and the boundary 
terms in the following way (^ is free parameter, irrelevant for the spectrum) 

A, = -l/2(gV2 + ,-1/2)^ h ^ i/2(,i/2 _ ^-1/2)^ 
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Bi = — — (a + 7 + (a - -f)af - 2afia^ - 2jfi ^af) , 

{(3 + 6-{P- 6)al - 25^9^/2-1/2^- _ 2/3^-1^-1^/2+1/2^+) 

For nonequilibrium systems, as opposed to the ones in equilibrium, the boundary conditions 
are of major importance. Emphasizing the dependence of the steady state behaviour on the 
boundary rates we study tridiagonal algebras of the simple exclusion processes which reveal 
deep algebraic properties of the latter. In a recent paper |14j we have shown that the boundary 
symmetry of the open partially asymmetric exclusion process is an Askey-Wilson algebra, the 
coideal subalgebra of Uq{s'u{2)), which allows for the exact solvability in the stationary state. 
From the boundary AW algebra the tridiagonal boundary algebra of the ASEP follows through 
the natural homomorphism. The particular values of the structure constants determine the bulk 
tridiagonal algebra of the ASEP. It is generated by the matrices Dq and Di of the matrix product 
ansatz and the defining relations have the form of the level zero Uq{su{2)) g'-Serre relations. 
This suggests a framework for analysis of the asymmetric exclusion process. Namely, given 
the Uq{su{2)) i?-matrix operator R{zi/z2) E Endcl4i Vz2, where Vz is the two-dimensional 
Uq{su{2)) evaluation module, satisfying the Yang-Baxter equation 

Rl2{zi/ Z2)Rl2.{zi)R-22.{z2) = R2?.{z2) Ri:i{zi) Ri-2{zi / Z2) , 

then the transition matrix F (or equivalently the Hamiltonian) is written as -ff = ^ Hu+i, where 
the two-site density is obtained as 



Hii+l — -j-PRii+i 
du 



n=0 



with P the permutation operator and 21/2:2 = e". The generators act on the quantum space by 
means of the infinite coproduct and the invariance with respect to the affine Uq{su{2)) manifests 
in the property 

[F,A~(G,.)]=0 

for any of the generators G^ of Uq{su{2)). If we introduce for finite chain a boundary of 
a particular form, such as diagonal boundary terms, the symmetry is reduced to Uq{su{2)) and 
the invariant Hamiltonian [13] -ff^xz i^ known [12] to describe a bulk diffusive system with 
reflecting boundaries. In the presence of a general boundary the symmetry is further reduced 
to a boundary symmetry whose generators are constructed as linear covariant objects with 
respect to the bulk quantum group. It turns out that the boundary symmetry of the asymmetric 
exclusion process is the AW algebra whose elements possess the coproduct properties of two-sided 
coideals of the bulk quantum Uq{su{2)) symmetry. We have implemented [TJ] the generators of 
the AW algebra to construct an operator valued i^-matrix, a solution to the spectral dependent 
boundary Yang-Baxter equation, also known as a reflection equation 

Rizi/z2){K{zi) <g) I)R{ziZ2){I <E) K{z2)) -{I® K{z2))R{ziZ2){K{zi) ® I)R{zi/z2) = 0. 

The K{z) matrix is the basic ingredient of the inverse scattering method [15] and such a connec- 
tion points out to exact solvability beyond the stationary state and description of the dynamics 
of the process. This can be achieved by applying the technique of the quantum inverse scattering 
method, and in particular the underlying Bethe Ansatz method to exactly calculate some quan- 
tities of physical interest, such as the common spectrum of the commuting conserved quantities. 
In the following we first review the algebra of the asymmetric exclusion process, both the 
boundary and the bulk one. The boundary AW algebra is naturally mapped to a tridiagonal 
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algebra whose irreducible infinite-dimensional modules are the AW polynomials as well. We 
study the q = 1 limit of the tridiagonal algebra and show that it appears to be the bound- 
ary symmetry of the symmetric exclusion process. The structure constants of the tridiagonal 
boundary algebra are given in terms of the parameters of the symmetric exclusion process. From 
the boundary algebra, analogously to the case of the ASEP, one obtains, for particular values 
of the structure constants, the bulk tridiagonal algebra of the symmetric process. We consider 
the Askey-Wilson algebra of the totally asymmetric exclusion process which can be viewed as 
a particular q = limit of the Askey-Wilson algebra. The consequences of these symmetry 
properties for the exact solvability of the simple exclusion model in the stationary state within 
the matrix product approach are discussed. 

2 The model within the matrix product approach 

The matrix product approach (MPA) was developed with the aim to describe the stationary 
behaviour of many-body systems with stochastic dynamics. The idea is that the steady state 
properties of the ASEP can be obtained exactly in terms of matrices obeying a quadratic algeb- 
ra |16lll7j . For a given configuration (si, S2, . . . , sl) the stationary probability is defined by the 
expectation value 

p/ N _ {w\Ds^Ds2---Dsl\v) 

where Dg. = Di, if a site i = 1,2,... ,L is occupied and -Dg. = Dq, if a site i is empty. 
The quantity Z^ = {w\{Dq + Di) \v) is the normalization factor to the stationary probability 
distribution. The operators Di, i = 0,1 satisfy the quadratic (bulk) algebra (known as diffusion 
algebra [IB]) 

DiDo - qDoDi = xiDq - DiXq, xq + xi = 

with boundary conditions of the form 

{PDi - dDo)\v) = xo\v), {w\{aDo - 7D1) = -{w\xi 

and (wlv) ^ 0. The boundary conditions define the two vectors (u;| and \v) which enter the 
expectations values for the stationary weights. The open model with boundary processes depends 
on five parameters. These are the bulk probability rate and the four boundary rates. The 
partially asymmetric exclusion process (PASEP) corresponds to < g < 1. The limit cases q = 1 
and q = are the symmetric exclusion process (SSEP) and the totally asymmetric exclusion 
process (TASEP). In the case of TASEP the boundary processes depend on two boundary rates 
corresponding to incoming particles at the left end of the chain and outgoing particles at the 
right end. 

Given the representations of the quadratic algebra and the boundary vectors, one can evaluate 
all the relevant physical quantities, such as the mean density at a site i, i = 1, . . . ,L, (sj) = 
'"'^^ " — — ^^^ " — —, correlation functions, the current J through a bond between site i 

and site i -|- 1 which has a very simple form J = — §^ . In most studied examples one uses infinite- 
dimensional representations of the quadratic algebra. Finite-dimensional representations |lll [T^ 
have been considered too and they simplify calculations. Due to a constraint on the model 
parameters they define an invariant subspace of the infinite matrices and give exact results only 
on some special curves of the phase diagram. 

Exact results for the ASEP with open boundaries were obtained within the MPA through 
the relation of the stationary state to g-Hermite [20] and Al-Salam-Chihara polynomials [21] in 
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the case 7 = (5 = and to the Askey-Wilson polynomials [22] in the general case. The MPA 
was readily generalized to many-species models [23] and to dynamical MPA |24| . 

Emphasizing the dependence of the open ASEP in the stationary state on the boundary rates 
and the equivalence with the spin 1/2 XX Z chain we can represent |25j the boundary operators 
in the form 

— 5 OL — T 

\-q \-q 

where yl, A* are linear combinations of the Uq{su{2)) generators in the evaluation representation 

^ = -VT=q' ^^-VT=q^-' T^^' ' 

The operators N, A± generate the Uq{su{2)) algebra {0 < q < 1) 

[N,A±] = ±A±, [A_,A^] = ^j^-^ (1) 

with a central element Q = A^A- — '^ . ^,^ ~'^-ii2\2 • The pair of operators A, A* satisfy the 
relations of the boundary Askey-Wilson algebra 

[[A, A%, A]q = -pA* -uA-rj, [A*, [A, A%], = -p*A - ujA* - 7?* (2) 

with structure constants 

p = xl(35q~\q'/' + ,-V2)^ p* = xla^q~Hq'/' + q^'/^ , (3) 

-iv = xliP -5){^-a)- xliP^ + a6){q'/^ - q-'/^)Q, 

V = ,'i\,'i' + ,-'^')4 {mi - -)Q + ^^;v?i%"'^ ' 

,* = qy\q^l^ + q~y^)xl (a,(/3 - 5)Q + ^" ;,,I^l;^.t f ^ 

Throughout the text [X, "Kjg = q^/'^XY — q~^''^YX. We note the relative difference of the 
factor q^'"^ in p (q'~^' ^ in p*) compared to [2^ which is simply due to a rescaling of the genera- 
tors A, A* . It is important to stress (see next section) that the Askey-Wilson algebra is defined 
up to an affine transformation of rescaling the generators by scalars t, t* . This property can 
be used to even hide the parameters (35 in p (07 in p*). It is known [261 [27] that the infinite- 
dimensional irreducible modules of the AW algebra are the AW polynomials depending on four 
parameters a, 6, c, d. In the basic representation one A* is represented by a diagonal and ^4 - by 
a tridiagonal matrix. The sign of p* (respectively of p in the dual representation) is essential [2T| 
for the spectrum of the diagonal operator. Since for a representation of the AW algebra we use 
the Uq{su{2)) algebra, which is the limit case of a (u, —u), u < algebra, the spectrum of the 
diagonal operator is of the form sinh (see [2^ for the details). The boundary conditions of the 
ASEP uniquely relate the four parameters of the AW polynomials to the four boundary rates 
a = K^(a,7), 6 = kj^{P,S), c = K*_{a,^), d = k-{P,S), where Kj! (z^, r) is 



) [u - T - {1 - q)) ± y/{i^ - T - {I - q)Y +"4^ 



2v 
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If we use the {u, u) algebra, the Hmit case of which corresponds to q + q~ in the nominator 
of equation ([1]), then the spectrum of the diagonal operator is of the form cosh and tv^ [u, r) — > 
KW(z.,-r). 

The transfer matrix Dq + Di and each of the boundary operators generate isomorphic AW 
algebras |25j . In the tridiagonal representation the transfer matrix Dq + Di satisfies the three- 
term recurrence relation of the AW polynomials, which was explored in [22] for the solution of 
the ASEP in the stationary state. The exact calculation of all the physical quantities, such as 
the current, correlation functions etc, in terms of the Askey-Wilson polynomials was achieved 
without any reference to the AW algebra. The exact solution in the stationary state in terms of 
the AW polynomials was, in our opinion, possible due to the appearance of the AW algebra as the 
boundary hidden symmetry of the ASEP with general boundary conditions. We have constructed 
an AW algebra operator valued iT-matrix [T3], a solution to the boundary Yang-Baxter equation. 
The relation of the boundary AW algebra to the JiT-matrix, satisfying the reflection equation, 
reveals deep algebraic properties of the ASEP allowing to extend the exact solvability beyond 
the stationary state. The boundary Yang-Baxter equation is the basic ingredient of the inverse 
scattering method [15] and one can use its solutions for description of the dynamics of the 
process. With the help of the i^(z)-matrix one can apply the technique of the quantum inverse 
scattering method to construct the spectral dependent transfer matrix and then implement the 
underlying Bethe Ansatz method to exactly calculate some quantities of physical interest, e.g. 
the exact spectrum of the transition matrix F, or equivalently the Hamiltonian. 

3 Bulk tridiagonal algebra of the simple exclusion process 

For completeness we first recall the bulk algebra of the partially asymmetric exclusion process 
(PASEP). 

The bulk tridiagonal algebra of PASEP: The operators Dq, Di of the partially asymmetric 
exclusion process and their (/-commutator [DQ,Di]q form a closed linear algebra 



[DQ,Dl]g=D2, 



V2^, ^nVa _ ^-1/2n 



pi, [Dq, D,],]g = q-'/'x^{q'/' - q-'/'){DQ, D^] 

- q-^xJDQ + q~^XQXiDi - XQq~^'^{q^'^ - q^^'^)Dl, 
[[DQ,Di]q,DQ]g = -XQq-^'\q^'^ - q-^'^){DQ,Di} 

- xlq-^Di + XQXiq-^DQ - x^q-^^q^/^ - q-'^'^)Dl, 

where {Dq,Di} = DqDi + DiDq. The algebra can equivalently be described as a two-relation 
algebra for the pair Dq, Di 

DqDI -{q + q~^)DiDQDi + DJDq + Xiq-^/^{q^'^ - q~^'^){DQ, Di} 

= xlq~^DQ - XQXiq~^Di + XQq~^'\q^'^ - q~^'^)Dl, 
DlDi -{q + q-^)DQDiDQ + D^DI - XQq-^'^{q^'^ - q~^'^){DQ, A} 

= xlq~^Di - XQXiq~^DQ + Xiq~^/^{q^'^ - q~^'^)Dl. 

The above relations are the well known Askey- Wilson relations 

A^A* -{q + q-^)AA*A + A*A^ - -^{AA* + A* A) = pA* + tM^ + w^ + 77, 

A*'^A -{q + q-^)A*AA* + AA*^ - j*{AA* + A* A) = p*A + jA*^ + uA* + r]*. (4) 

The algebra (JH) was first considered in the works of Zhedanov |271 [28] who showed that the 
Askey-Wilson polynomials give pairs of infinite-dimensional matrices satisfying the Askey- 
Wilson (AW) relations. It is recently discussed in a more general framework of a tridiagonal 
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algebra [291 EO] , that is an associative algebra with a unit generated by a (tridiagonal) pair of 
operators A, A* and defining relations, obtained by taking the commutator with A, A* in the 
first (second) line of ^ respectively 

[A, A^A* - (3AA*A + A*A^ - -^{AA* + A* A) - pA*] = 0, 
[A\A*'^A - pA*AA* + AA*^ - j*{AA* + A* A) - p*A] = 0. 

In the general case a tridiagonal pair is determined by the sequence of scalars (3, 7, 7*, p, p* . 
(We keep the conventional notations, used in the literature, for the scalars of a tridiagonal pair; 
P and 7 should not be confused with the ASEP boundary rates.) Tridiagonal pairs have been 
classified according to the dependence on the scalars [29]. The example which is important for 
the present study are the Dolan-Grady relations [31] with P = 2, ^ = j* = 0, p = k"^ , p* = k*'^ 

[A,[A,[A,A*]]]=k^[A,A*], [A*,[A*,[A*,A]]] = k*^[A*,A]. 

Tridiagonal pairs are determined up to an affine transformation 

A^tA + c, A* ^t*A* + c*, 

where t, t* , c, c* are some scalars. The affine transformation can be used to bring a tridiagonal 
pair in a reduced form with 7 = 7* = 0. 

As seen from the Askey-Wilson relations ([4]) for the matrices Dq and Di we have explicitly 

2-1 * 2 -1 -1 

p = XiQ , p =XQq , LU = -XQXiq , 

-1/2/ 1/2 -l/2\ * -1/2/ 1/2 -l/2\ 

J = -xiq ^i\q''''-q ''^), 7 = x^q ''\q^i''-q '-''') 

and r] = rj* = 0. Besides 7 = 7*, p = p* due to xq + xi = 0. We can now rescale the 
operators Dq, Di to set 7 = 7* = 0. This is achieved with the help of the transformations 



ql/2_q-l/2' ^ ^ ql/2_q-l/2- 

However the shift of the generators amounts to a tridiagonal pair with sequence of scalars 
P = —{q + q~^), 7 = 7* = 0, p = /9* = 0. Thus the operators of the ASEP matrix product 
ansatz obey the relations of a tridiagonal algebra 

[Di,DoDl -{q + q-^)DiDoDi + D^Do] = 0, 

[Do, D^Dl -{q + q-^)DoD^Do + D^D^] = 0, (5) 

which is a special case of the tridiagonal relations of the ASEP boundary operators with P = 
q + q~^ and 7 = 7* = /) = p* = 0. These relations are the g-Serre relations for the level 
zero Uq{su{2)) adjoint representation. As already mentioned we have defined the Askey-Wilson 
algebra with two generators through a homomorphism to the quantized affine Uq{su{2)) [T3j. It 
is important to understand the role of the parameters xq, xi in the quadratic algebra of the open 
ASEP from the point of view of this homomorphism. These parameters can be interpreted as 
the weights associated to the Cartan generators i/j, i = 0, 1 of the Uq{su{2)) finite-dimensional 
module with weight space W^ = (z/ G W\q^^v = q^^u). The Uq{s'u{2)) module W has level k 
if the central element q^o+Hi ^^^^g qj^ ^^ ^g ^^iq scalar g'^ [33]. Hence, it follows from the bulk 
tridiagonal algebra dSj) that the matrices Dq, Di of the quadratic bulk algebra of the open ASEP 
with xq + xi correspond to the two-dimensional level Uq{s'u{2)) adjoint module. 

The bulk Askey-Wilson algebra of the symmetric exclusion process follows immediately as 
the limit (7 — > 1 of the bulk tridiagonal algebra of the asymmetric exclusion process. Hence we 
have: 
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Proposition 1. The operators Dq, Di of the symmetric simple exclusion process and their 
commutator [Dq,Di\ form a closed linear algebra 

[Do, Di] = D2, [Di, [Do, Di]] = xi{Do, Di} - xJDo + xox^Di - xoDf, 

[[Do, Di],Do] = -xo{Do, Di} - x^Di + xoXiDq - xiD^. 

The proposition can be independently verified by directly using the explicit form of the MPA 
quadratic relation for the case of symmetric diffusion. The algebra can equivalently be described 
as a two-relation algebra for the pair Do, Di 

DoDf - 2DiDoDi + DJDo + xi{Do, Di} = xJDo - xoXiDi + xoDJ, 
DlDi - 2D0D1D0 + DiDi - xo{Do, Di} = xIDi - xoXiDo + x^dI 

As seen from the explicit form of the algebra, the matrices Do, Di of the symmetric exclusion 
process satisfy the Askey- Wilson relations with the sequence of scalars 

/3 = 2, P = x1, P* = xl, io = -xoxi, 7 = -xi, 7* = -xq. 

We can transform the generators Do, Di by the affine shifts 

Do^ Do + xo, Di ^ Di- XI. 

By a straightforward calculation we again observe the same property as in the case of the 
partially asymmetric process. Namely, the shifts of the generators amount to a tridiagonal pair 
with the sequence of scalars 7 = 7* = 0, p = p* = 0. Thus the operators of the symmetric 
matrix product ansatz obey the relations of a tridiagonal algebra 

[Di,DoDl - 2D1D0D1 + DlDo] = 0, [Do, D^dI - 2DoDiDo + dId^] = 0, 

which are the Dolan-Grady relations with A; = /c* = 0. In the next section we show that the bulk 
tridiagonal relations are the special case of the tridiagonal relations obeyed by the boundary 
operators of the symmetric process. 

4 Boundary algebra of the symmetric exclusion process 

We consider the symmetric simple exclusion process (SSEP) with most general boundary con- 
ditions of incoming and outgoing particle at both ends of the chain. Within the matrix product 
ansatz the quadratic bulk algebra of the SSEP is the q = 1 limit of the deformed quadratic 
algebra of the ASEP. In the previous section we derived the bulk tridiagonal algebra of the 
symmetric process and it turned to be the g = 1 limit of the bulk tridiagonal algebra of the 
asymmetric process. 

We are going to obtain the boundary algebra of the symmetric process as the g = 1 limit of 
the tridiagonal boundary algebra of the asymmetric simple exclusion process. We start with the 
observation about the natural homomorphism of the tridiagonal algebra (TA) generated by the 
pair A, A* and the Askey- Wilson algebra (AW) of the ASEP defined by equation ([2]), namely 
TA — > AW . As already mentioned, this is readily verified by taking the commutator with A, 
respectively A* of the first line, respectively the second line of equation ([2]), which gives 

[A,[A,[A,A%,\^-,]]=p[A,A*], [A\[A*[A\A\]^-A=p*[A*,A] (6) 

with p, p* depending on the five parameters of the ASEP, as given by ^. The operators A, A* 
were introduced as the shifted boundary operators of the ASEP. The limit q' = 1 of the tridiagonal 
algebra ^ provides a way to determine the boundary algebra of the SSEP. 
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The boundary operators of the SSEP can be represented in the form (i.e. hnear combinations 
in terms of the level zero affine su{2) generators) 

(3Di — 5Dq = —xij3Aj^ — xq5A^ — {xiP + xo6)N — xifi — xq5, 

uDq — ^Di = xquAj^ + xijA- + {xoa + xi^)N + xqo + xij, (7) 

where A±, N are corresponding operators in the limit g — > 1 of equation ([H). We separate the 
shift parts from the boundary operators. Denoting the corresponding rest operator parts by A 
and A* we write the left and right boundary operators in the form 

f3Di — 5Dq = a — xi(3 — xq5, oDq — ^Di = A* + xqu + 0:17. (8) 

Then we have the following 

Proposition 2. The operators A and A* defined by the corresponding shifts of the boundary 
operators of the open symmetric exclusion process 

A = f3Di - 6D0 + {xiP + xo5), A* = aDo - 7D1 - (xoa + 2:17) 
and their commutator 

[A, A*] = AA* - A* A 
form a closed linear algebra, the boundary tridiagonal algebra of the SSEP 

[A, [A, [A, A*]]] = p[A,A*], [A*, [A*, [A*, A]]] = p*[A*,A], (9) 

where the structure constants are given by 

—p = xqXi(3S, —p* = XQXia^ (10) 

This proposition is straightforward to verify by taking the g — > 1 limit in the chain of 
homomorphisms TD -^ AW -^ Uq{su{2)) or independently by directly using equation ([7]). 

As readily seen from the definition ([8]) the (shifted) boundary operators of the symmetric 
exclusion process obeying the algebra ([9]) form a tridiagonal pair with /? = 2, 7 = 7* = 0, 
and p, p* given by ()10p . The tridiagonal boundary algebra of the symmetric process is the limit 
(7 = 1 of the deformed boundary algebra of the ASEP as the irreducible modules of the algebra 
in the symmetric model, i.e. the Wilson polynomials, are the q = 1 limit of the Askey- Wilson 
polynomials. The important properties of the deformed Askey-Wilson algebra remain valid in 
the proper limit g — > 1. Stated more precisely which will correspond to the historical development 
of these algebras, for generic q the deformed Askey-Wilson algebra is the g-generalization of 
the Onsager algebra in the form of the Dolan-Grady relations. Its irreducible modules, i.e. the 
Askey-Wilson polynomials are g'-counterpart of the Wilson polynomials. For applications, howe- 
ver one can even use the limit cases of the Wilson polynomials. These are the continuous Hahn 
and dual Hahn polynomials (see [M] and [35] for details on these polynomials). In both cases 
there exist limiting procedures to further obtain the Meixner-Pollaczek polynomials from which 
the Laguerre polynomials can be obtained. Let Pn{x,(f>) denote the n-th Meixner-Pollaczek 

The Laguerre polynomials can be obtained by the substitution p = 1/2(A + 1), x — > —l/2(j)~^{x) 
and letting (/> — > 



JL-o^«'""'^(-2^^*)=^i-'(^) 
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By definition the Laguerre polynomials have the form 



n! "^ V A + 1 
with orthogonality condition 

r e-^xHlLidx = Ii!i±A±i)5_, A > -1 

Jo n\ 

and recurrence relation 

(n + l)Li2,(x) - (2n + A + 1 - x)lW(x) + {n + X)l[%{x) = 0, 

where L_{(x) = and Lq (x) = 1. One can identify 

a + /? + 7 + J 
"(a + 7)(/3 + <5) ■ 

Denoting 

we rewrite the orthogonality condition in the form 

1 f-OO 

The vectors \l{x)) = (/o(a^),^i(0), . . . )* and {l{x)\ = {lo{x),li{0), . . .) form the basis for the 
tridiagonal and the diagonal representation of the generators (and for the dual one). As it was 
proved to be the case for the ASEP [52], each boundary operator of the symmetric process, 
together with the transfer matrix operator Dq + Di, forms an isomorphic AW algebra whose 
tridiagonal representation follows from the three-term recurrence relation, consistent with the 
orthogonality condition 

{Do + Di)\l{x))=x\l{x)), {l{x)\{Do + Di) = {l{x)\x. (11) 

These properties can be used to exactly calculate the physical quantities [22] in the stationary 
state. The results [2] for the partition function Zi and the current are reconstructed 



r(A + l) ' X + L 

The one point function [2] 

a 1 aB — 7(5/1 \ 

in) = —— - .^r I / MR^s\ ( —— + ^ - 1 ' i = l,2,...,L 

a + 7 A + L (a + 7)(/^ + 0) va + 7 / 

shows that the particle density has a linear profile. Without reference to the AW algebra the 
eigenvalue equation (llip was used in [22] to calculate the physical quantities of the symmetric 
exclusion process, as known from the matrix product approach. In our opinion it is the boundary 
tridiagonal algebra of the symmetric exclusion process that allows for the exact solvability of the 
symmetric process in the stationary state and leads to a generalization of the matrix product 
method. 
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5 Nonlocal conserved charges 

of the symmetric exclusion process 

In the previous section we have shown that the boundary symmetry of the symmetric exclusion 
process is the tridiagonal algebra ([9]) with the sequence of scalars /3 = 2, p, p* . This algebra is 
the q = 1 limit of the tridiagonal algebra mapped through the natural homomorphism to the 
Zhedanov algebra 74Ty(3) defined for < g < 1. The exact calculation of the physical quanti- 
ties in the case of the symmetric exclusion process has been obtained in terms of the Laguerre 
polynomials implementing [22] the three-term recurrence relation to define the tridiagonal rep- 
resentation of the transfer matrix Dq + Di. As we pointed out, the ultimate relation of the 
exact solution in the stationary state to the Laguerre polynomials was possible due to the q = 1 
boundary hidden symmetry of the SSEP with general boundary conditions. 

The defining relations ([9]) of the q = 1 tridiagonal boundary algebra can (with p, p* expressed 
in terms of the boundary rates according to (jlOp ) are the well-known Dolan-Grady relations 
for the shifted boundary operators. The importance of the Dolan-Grady relations is that given 
a self-dual Hamiltonian 

H = fA + f*A*, 

where /, /* are some coupling constants, and A, A* satisfy the relations ([9]), then one can 
construct (an infinite set of) conserved commuting charges 

Q2n = f{R2n — R2n~2) + f*{R2n — R2n~2) 

in terms of the quantities 

R2n = --[^[A*,R2n-2]] " ^2„-2 (12) 

P 

with Rq = a, Qq = H. In the case of the symmetric exclusion process from the boundary 
conditions we have 

{w\A — {xif3 + xo5)\v) = xo{'w\v), {w\A* + {xqo + xi7)|t)) = —xi{'w\v. 
Hence 

{w\xo{A* + {xoa + X17)) +xi{A- (xo/3 + xi6))\v) = 
and we can interpret the quantity 

xo{A* + {xoa + X17)) +xi{A- {xoj3 + xiS)) 

as the Hamiltonian of the symmetric simple exclusion process in the auxiliary space. As we 
know the generators of the tridiagonal algebra are determined up to shift transformations. In 
view of this property it seems more convenient to consider the shifted Hamiltonian 

Hs = xqA* + xiA, 

which is self-dual, if we define x^ = xq. Then we can straightforward apply the prescription of 
Dolan and Grady to obtain the conserved nonlocal charges. Taking into account the shifts we 
find the result 

Rq = A — xiP — xo5, Rq = A* + Xoa + X17, 
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2 

R2 = — [A — xif3 — X06, [A* + xoa + xi^,A — xiP — xo5]] — A* + x^a + X17 
P 

and so on, according to formula (J12p . The expressions A — xi[3 — xq5 and A* + xqu + 3:17 are 
the right and left boundary operators 

B^ = /3Di-6Do, B^ = aDi--fDo 

respectively which acquire a very important physical meaning. The boundary operators satis- 
fying 

[B«, [B^, [B^,B^]]] = -xoXipd[B^,B^], 

[B\ [B\ [B\ B%] = -xoXiaj[B\ B^] (13) 

are the nonlocal conserved charges of the symmetric exclusion process with the help of which 
the (infinite) set of conserved quantities of the process are constructed. By properly rescaling 
the operators one achieves equal coefficient factors on the RHS of equations (1130 . We note that 
quantum integrals of motion for the XXX Heisenberg infinite chain (known to be related to 
the symmetric exclusion process) were first obtained in [32]. For the open symmetric exclusion 
process the existence of nonlocal conserved quantities will result in the exact solvability of the 
system beyond the stationary state due to the boundary hidden tridiagonal symmetry of the 
symmetric exclusion process. 

6 Askey— Wilson algebra of the totally asymmetric 
exclusion process 

The tridiagonal and the AW algebra of the totally asymmetric exclusion process cannot be 
obtained directly as the limit q = of the partially asymmetric process. The procedure is more 
involved. We derive this algebra from the quadratic algebra of the totally asymmetric process. 
We start with the quadratic algebra 

DiDo = Di + Do. (14) 

From this algebra the following relations follow 

DiDoDi = Dl + DoDi, DqDiDo = DoDi + D^ (15) 

and 

DlDo = Dl + Di + Do, DiD^ = Di + Do + Dl (16) 

which can alternatively be written as 

DiDoDi - DlDo = [Do, Di], D^D^Do - D^dI = [Do, Di]. (17) 

The LHS of (jl7|) are respectively 

[DiDo,Di], [Do,DiDo]. (18) 

Hence we have 

D^[Do,Di] = [Do,Di], [Do,Di]Do = [Do,Di]. (19) 
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Examples of matrices obeying the above relations are given by the equation ^ and (36), (38) 
in [12] ■ We can shift the operators Dq, Di by 1 (or respectively by the constants cq = a, ci = b) 

Do^Do + 1, Di^D, + 1. 

Then we can write subsequently DiDq = Z, where either Z = 1 or Z = ab. (In the case of 
the totally asymmetric exclusion process within the matrix product approach ab = a(3, where a 
and (3 are the probability rates for the particles to be added and removed at both sides of the 
linear chain.) Hence 

DiDoDi = ZDi, DoDiDo = ZDo, (20) 

DfDo = ZDi, DiDl = ZDq. (21) 

Relations (I17p for the shifted generators become 

D^DoDi - DlDo = 0, DoD^Do - D^D^ = (22) 

and consequently 

Di[Do,Di]=0, [Do,Di]Do = 0. (23) 

We need to emphasize that multiplying in equations ()22p . the first equation by Di from the 
right and the second equation by Dq from the left, we obtain 

DiDoDJ - DlDoDi = 0, DIDiDq - D^DiDl = 0, (24) 

which define the q = Q limit of the g-Serre relations of Uq{sl{2)), i.e. the q = limit of the level 
zero Uq{sl{2)) adjoint representation. This is consistent with the definition of the Askey-Wilson 
algebra such, that it yields a deformation of the level zero Uq{sl{2)) g-Serre relations (see [H] 
for details). However, we have now the additional relation 

It is important to emphasize that the matrices obeying (fT^ - (fT9]) and (f20]) - (i23]) are upper bi- 
diagonal and lower bidiagonal. It is known that the basic representation of the AW algebra can 
be equivalently considered as a representation on the space of c-number sequences where one 
of the generators acts as a tridiagonal operator, while the other generator acts as a diagonal 
operator. In our case the analogues of the diagonal and tridiagonal matrices in the g = limit 
of the basic representation are given by [Dq, Di] and Di + Dq respectively. 

Proposition 3. The q = ^ AW algebra, as following from the quadratic algebra ()14p . depending 
on only two constants a, b (i.e. ei = a + b, 62 = ab) in the basic representation, is defined by 

DiDQ = e2, Di[Dq,Di]=0, [Dq,Di]Dq = 0, 

[DQ,Di]iDQ + Di)[DQ,Di] = 0, i[DQ,Di])^ = e2[DQ,Di]. 

The matrix Dq+Di, denoted hereafter D, plays the role of the tridiagonal operator, and [Dq,Di], 
denoted hereafter D* , - the role of the diagonal operator, in the basic representation. 
Thus we have 

[Dq, Di]{Dq + Di)^[Dq, Di] = el[DQ, D^], 

[Dq,Di]{Dq + Di)[Dq,Di] = 0, {[DQ,Di])^ = e2[DQ,Di]. 

We can now shift the matrix {Dq + Di) = D^D + a + b to obtain a tridiagonal matrix with 
entries on the main diagonal too. 
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Definition 1. The q = limit of the AW algebra depending on only two constants a, b, where 
ei = a + b, 62 = ab, is generated by a tridiagonal operator D and a diagonal operator D* with 
defining relations in the basic representation 

D*DD* = eie2D*, {D*f = e2D*. (25) 

One can alternatively consider defining relations for g = limit of the AW algebra in a rep- 
resentation associated with the g = limit of the level zero adjoint Uq{sl{2)) (equation ()24p ). 

Definition 2. The q = Q limit of the AW algebra, generated by upper diagonal and lower 
diagonal matrices Di and Dq and depending on only two constants a, 6, is defined by 

DiDoDi=abDi, DoDiDo = abDo. (26) 

The boundary tridiagonal AW algebra of the totally asymmetric exclusion process corresponds 
to the values a = a, b = [3 while the bulk AW algebra is obtained for a = b = 1. 

We note that the two representation dependent definitions of the algebras considered above 
have in common the first relation in formula (j25p and both relations in (j26p . which can be 
unified as 

ABA = c{a,b)A, (27) 

where A, B are the generators of these algebras and c{a, b) is a constant depending on the 
parameters a, b. We can multiply equation (j27p by B subsequently from the left and from the 
right to obtain 

ABAB - BAB A = (28) 

if c{a, 6) = or 

ABAB ~ BAB A = c{a, b) {AB - BA) (29) 

if c[a^ b) 7^ 0. We can consider equation (j29p as defining an alternative g = limit algebra, 
obtained through the natural homomorphism to the g = limit Askey-Wilson algebra, with 
structure constant c{a, b). These relations can be very useful for applications. Namely, with the 
additional condition A^ = A, which is a Hecke type relation equations (j28p and (j29p have the 
form of a reflection equation and a modified boundary Yang-Mills equation, respectively. The 
interpretation of equations (|28p and (j29p as reflection equations should be associated with the 
proper iZ-matrix operator depending on a parameter t ^ q. 

Thus for the totally asymmetric exclusion process we obtain the bulk algebra 

DiDoDiDo - DoDiDoDi = DiDq - D^Di 
and the boundary algebra generated by the right D^ and left D*L boundary operators 

D^ = l3Di, D^ = aDo 
subject to the relations 

D^D^D^D^ - D^D^D^D^ = a(3{D^D^ - D^D^). (30) 

It is worth studying the connection of the exact solvability of the totally asymmetric exclusion 
process in the stationary state to the integrability properties based on the boundary Yang-Baxter 
equation. The consequences of equation (j30p will be to naturally extend the exact solvability 
beyond the stationary state. In [36] Bethe ansatz equations were derived and the exact spectrum 
of the transfer matrix of the totally asymmetric exclusion process was analyzed. In our opinion 
the boundary AW algebra is the hidden symmetry behind the Bethe ansatz solvability of the 
totally asymmetric exclusion process. 
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7 Conclusion 

We have studied the symmetry properties of the simple exclusion process. Our consideration ex- 
tends the previously obtained results about the AW algebra (known as the AW{3,q)) Zhedanov 
algebra) to be the boundary hidden symmetry of the asymmetric exclusion process. We have 
shown that the boundary symmetry algebra of the symmetric simple exclusion process is the 
q = 1 limit of the boundary tridiagonal algebra of the partially asymmetric exclusion process. It 
is this algebra, to be denoted TA{q = 1), that allows for the exact solvability of the symmetric 
model in the stationary state. The defining relations of the algebra TA{q = 1) are the known 
Dolan-Grady relations. The consequence of this fact is that one can define conserved nonlocal 
charges of the symmetric process, which allows for extending the exact solvability beyond the 
stationary state. We have also derived the bulk and boundary AW algebra of the totally asym- 
metric exclusion process, which can be viewed as a particular case of the algebra ^VF(3, q = 0). 
The defining relation of the boundary algebra of the totally asymmetric process leads to a 
boundary Yang-Baxter equation which might be interpreted again as the deep algebraic reason 
behind the exact solvability. 
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